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ABSTRACT 

This  proposal  requests  funding  for  acquisition  of  a  High  Performance  Computer  (HPC)  cluster  to  expand  and  strengthen  the  computational 
facility  of  W.  M.  Keck  Computational  Materials  Theory  Center  (CMTC)  at  California  State  University  Northridge  (CSUN).  The  proposed 
HPC  cluster  tremendously  enhanced  the  computational  capability  of  CMTC  and  allowed  the  researchers  to  tackle  challenging  materials 
problems.  The  HPC  cluster  has  been  and  will  be  used  to  perform  first-principles  simulations  to  predict  charge  carrier  mobilities,  exciton 
diffusion  and  interfacial  charge  separation  in  all-organic  and  hybrid  organic-inorganic  solar  cells.  The  outcome  of  the  project  is  the 
development  of  a  computational  capability  that  can  potentially  shift  the  paradigm  of  materials  design  based  on  time-consuming  trial-and- 
error  experiments  and  significantly  reduce  the  time  and  labour  required  for  materials  development.  The  proposed  cluster  will  also  play  an 
important  role  for  education  and  training  of  STEM  students  at  CSUN.  Not  only  can  the  students  involved  in  the  DoD  projects  use  the  cluster, 
but  the  students  taking  the  computational  materials  courses  can  also  access  the  cluster  for  their  coursework.  In  addition,  the  local  high  school 
teachers  and  students  attending  the  summer  camps  on  Computational  Materials  Research  at  CSUN  will  have  the  opportunity  to  use  the 
cluster  as  well. 
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Technical  Objectives 

This  proposal  requests  funding  for  acquisition  of  a  High  Performance  Computer  (HPC)  cluster  to  expand 
and  strengthen  the  computational  facility  of  W.  M.  Keck  Computational  Materials  Theory  Center 
(CMTC)  at  California  State  University  Northridge  (CSUN).  The  proposed  HPC  cluster  will  tremendously 
enhance  the  computational  capability  of  CMTC  and  allow  the  researchers  to  tackle  challenging  materials 
problems  that  are  presently  beyond  the  reach.  The  HPC  cluster  will  be  used  to  perform  first-principles 
simulations  to  predict  charge  carrier  mobilities,  exciton  diffusion  and  interfacial  charge  separation  in  all- 
organic  and  hybrid  organic-inorganic  solar  cells.  The  first-principles  simulations  validated  through 
experiments  could  accelerate  the  discovery  of  novel  organic  semiconductors  for  applications  in 
photovoltaics,  light- emitting  diodes,  flexible  field-effect  transistors,  lighting,  and  displays.  The  proposed 
research  involves  advancement  of  first-principles  simulation  methodologies,  development  of  user-friendly 
software,  validation  of  the  software  through  experiments,  gaining  fundamental  understanding  of 
underlying  processes,  translation  of  the  understanding  to  materials  design  rules,  and  providing  theoretical 
guidance  and  computational  screening  for  rational  design  of  materials.  The  outcome  of  the  project  is  the 
development  of  a  computational  capability  that  can  potentially  shift  the  paradigm  of  materials  design 
based  on  time-consuming  trial-and-error  experiments  and  significantly  reduce  the  time  and  labour 
required  for  materials  development. 

The  proposed  cluster  will  also  play  an  important  role  for  education  and  training  of  STEM  students  at 
CSUN.  Not  only  can  the  students  involved  in  the  DoD  projects  use  the  cluster,  but  the  students  taking  the 
computational  materials  courses  can  also  access  the  cluster  for  their  coursework.  In  addition,  the  local 
high  school  teachers  and  students  attending  the  summer  camps  on  Computational  Materials  Research  at 
CSUN  will  have  the  opportunity  to  use  the  cluster  as  well. 


Technical  Approach 

The  PI  has  developed  a  first-principles  method  that  can  predict  charge  carrier  mobility  in  organic 
semiconductors  as  a  function  of  chemical  structures,  presence  of  defects,  temperature,  electric  field,  and 
carrier  concentration  entirely  from  first  principles;  i.e.,  there  is  no  empirical  input  or  adjustable  parameter 
in  the  simulations.  More  specifically,  this  approach  is  based  on  the  time-domain  ab  initio  non-adiabatic 
molecular  dynamics  for  simulating  phonon-assisted  electron  transitions.  The  electronic  energy  levels  and 
their  transition  rates  are  determined  from  the  density  functional  theory  (DFT),  taking  into  account  of  both 
inter-molecule  and  intra-molecule  contributions  and  treating  dynamic  and  static  disorder  at  an  equal 
footing.  This  multiscale  framework  incorporates  quantum  mechanical  determination  of  electron  transition 
rates,  mesoscale  master  equations  for  carrier  mobility  and  continuum  models  for  current-voltage 
characteristics  tailored  for  experimental  measurements.  For  interfacial  exciton  dissociation,  we  have 
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developed  a  first-principle  approach  based  on  the  time-dependent  density  functional  theory  (TDDFT)  to 
describe  exciton  states,  including  energy  levels  and  many-body  wave  functions.  The  non-adiabatic  ab 
initio  molecular  dynamics  (MD)  is  used  to  determine  the  phonon-assisted  transition  rates  between 
localized  exciton  states.  In  conjunction  with  the  Monte  Carlo  and  surface-hopping  methods,  this  approach 
can  simulate  exciton  dynamics  in  organic  semiconductors  at  different  temperatures. 


Progress: 


1.  Charge  Separation  and  Exciton  Dynamics  at  Polymer/ZnO  Interface  from  First-Principles 
Simulations  (Wu  et  al.,  J.  Phys.  Chem.  Lett.  5,  2649,  2014) 

Charge  separation  and  exciton  dynamics  play  a  crucial  role  in  determining  the  performance  of  excitonic 
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(d)  Hot  CT,  3.4  eV  (e)  CBM,  0.0  eV 
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photovoltaics.  Using  time-dependent  density  functional 
theory  with  a  range-separated  exchange-correlation 
functional  as  well  as  non-adiabatic  ab  initio  molecular 
dynamics,  we  have  studied  the  formation  and  dynamics  of 
charge-transfer  (CT)  excitons  at  polymer/ZnO  interface. 
The  interfacial  atomic  structure,  exciton  density  of  states 
(DOS)  and  conversions  between  exciton  species  are 
examined  from  first-principles.  The  interfacial  atomic 
structure  and  exciton  DOS  provide  the  basis  to  understand 
the  interfacial  charge  separation.  The  conversions  among 
the  intermolecular,  intramolecular,  and  CT  excitons  are 
examined  in  detail.  The  exciton  dynamics  exhibits  both 
adiabatic  and  non-adiabatic  characters.  While  the 
adiabatic  transitions  are  facilitated  by  C=C  vibrations 
along  the  polymer  (P3HT)  backbone,  the  non-adiabatic 
transitions  are  realized  by  exciton  hopping  between  the 
excited  states.  We  find  that  the  localized  ZnO  surface 
states  lead  to  localized  low-energy  CT  states  and  poor 
charge  separation.  In  contrast,  the  surface  states  of 
crystalline  C60  are  indistinguishable  from  the  bulk  states, 
resulting  in  delocalized  CT  states  and  efficient  charge 
separation  in  polymer/fullerene  (P3HT/PCBM) 
heterojunctions.  The  hot  CT  states  are  found  to  cool  down 
in  an  ultrafast  time  scale  and  may  not  play  a  major  role  in 
charge  separation  of  P3HT/ZnO.  Finally  we  suggest  that 
the  dimensions  of  nanostructured  acceptors  can  be  tuned  to  obtain  both  efficient  charge  separation  and 
high  open  circuit  voltages. 


Figure  1:  Charge  density  of  CT1  (a),  CT2  (b),  CT3  (c), 
and  the  hot  CT  state  with  the  hole  on  Pi  (d).  The  red  and 
blue  color  represents  the  charge  density  of  the  hole  and 
the  electron,  respectively. 


2.  Electric  Field  Control  and  Effect  of  Pd  Capping  on  Magnetocrystalline  Anisotropy  in  Fe-Pd  Thin 
Films:  A  First-Principles  Study  (Ong,  et  al.,  Phys.  Rev.  B  89,  094422,  2014) 

Using  ab  initio  electronic  structure  calculations,  we  have  investigated  the  effect  of  electric  field  and  of 
heavy  metal  cap  of  Pd  on  the  magnetocrystalline  anisotropy  (MCA)  of  Fe-Pd  ultrathin  film.  Analysis  of 
the  energy-  and  k-resolved  distribution  of  the  orbital-character  of  the  minority  spin  band  reveals  that  the 
perpendicular  MCA  of  the  uncapped  film  mainly  arises  from  the  spin-orbit  coupling  (SOC)  between 
unoccupied  Fe  dxy  and  occupied  Fe  dx2-}’2  states.  On  the  other  hand,  the  SOC  between  the  Pd-  and  Fe- 
derived  cZ-states  yields  negative  contributions  to  the  MCA.  We  find  that  the  sensitivity  of  the  surface 
anisotropy  energy  to  the  applied  electric  field  is  18  fJ/(Vm)  and  is  due  to  changes  in  the  occupation  of  the 
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surface  Fe  atoms  dxr-yi  and  (dxz;  dyz)  orbitals.  We  demonstrate  that  the 
thickness  of  the  Pd  cap  has  a  dramatic  effect  on  the  MCA  and  can  even 
switch  the  magnetization  from  out-of  to  in-plane  orientation.  The 
underlying  origin  is  the  change  of  the  position  and  orbital  character  of  the 
spin-polarized  quantum  well  states  induced  in  the  Pd  cap  by  varying  its 
thickness.  These  results  have  important  implications  on  exploiting  heavy 
metals  with  large  spin-orbit  coupling  (Ru,  Pd,  Ta,  Pt,  or  Au)  as  contacts 
with  ferromagnetic  thin  films  to  tailor  the  magnetic  switching  of  spintronic 
devices  by  tuning  the  cap  thickness. 


Figure  2:  Energy-  and  k-resolved  distribution  of  the  orbital  character  of  the  minority-spin  band  of  the  9  ML  FePd  film. 
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